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We calculate the optical conductivity using several models for unparticle or scale-invariant mat¬ 
ter. Within a Gaussian action for unparticles that is gauged with Wilson lines, we find that the 
conductivity computed from the Kubo formalism with vertex corrections yields no non-trivial de¬ 
viation from the free-theory result. This result obtains because at the Gaussian level, unparticles 
are just a superposition of particle fields and hence any transport property must be consistent with 
free theory. Beyond the Gaussian approach, we adopt the continuous mass formulation of unpar¬ 
ticles and calculate the Drude conductivity directly. We show that unparticles in this context can 
be tailored to yield an algebraic conductivity that scales as with the associated phase angle 

between the imaginary and real parts of arctan ^ = 60° as is seen in the cuprates. Given the 
recent results[l-3] that gravitational crystals lack a power-law optical conductivity, this constitutes 
the first consistent account of the conductivity and the phase angle seen in optimally doped 

cuprates. Our results indicate that at each frequency in the scaling regime, excitations on all energy 
scales contribute. Hence, incoherence is at the heart of the power-law in the optical conductivity in 
strongly correlated systems such as the cuprates. 
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I. INTRODUCTION 

Among the signatures of non-Fermi liquid behavior, the fractional power law observed in the optical conductivity 
of optimally doped cuprates, (t{(jj) oc 5] with a = 2/3, stands as a key failure of the standard theory of 

metals. This power law is ubiquitous, having been observed first in YBa2Cu306+5 [6, 7] with a = 0.5 — 0.7 and more 
extensively in Bi2Sr2CaCu208+5 [4, 8] in the mid-infrared frequency range (w = 500 cm~^ to 5000 cm~^). Consistent 
with this power law is the observed phase angle of 7rQ;/2[4]. Realizing that some sort of scale invariant entities must 
be relevant, Anderson[9] computed the mid-infrared power law in the cuprates using Luttinger liquid Green functions. 
He argued that in the holon nondrag regime, the vertex term can be dropped out of the conductivity calculation, 
^ ^ f dx f dtG^{x,t)G^{x,t)e^‘^*, where G® and G^ are the electron and hole Green functions of a 2D Luttinger 
liquid, respectively. Although no specific calculation was provided, he posited that the Green function of 2D Luttinger 
liquids (an unknown state of matter) is identical to its ID counterpart. Using the scaling properties of the product 
of Green functions, Anderson found the optical conductivity to have a power law of a form, a{uj) oc where 

a is the Luttinger parameter. 

Given that power laws are the finger print of criticality, Anderson’s starting point of scale-invariant excitations [9] is 
certainly reasonable. However since there is no basis for using ID Luttinger liquids in 2D, it is advantageous to resort 
to a relatively new tool which is capable of describing quantum critical matter in higher dimensions. Previously, one 
of us adopted[10] the unparticle picture of Georgi’s[ll] as a model for treating the breakdown of the particle picture 
in the normal state of the cuprates. The key idea here is that scale invariance places severe restrictions on the form of 
the propagator that are completely dimension independent and hence obviates the need to tether the description to 
ID Luttinger-liquid theory. We investigate two scale invariant models here. In the first model, we use a Gaussian form 
for the unparticle action as a basis for calculating the optical conductivity. We apply the vertex couplings of scalar 
unparticle and gauge fields introduced in Refs. [12, 13] and linear response theory at finite temperature. We show 
that the conductivity of unparticles is a constant factor multiplied by the conductivity of free massless scalar fields. 
In the second model, we turn to the continuous mass formalism that one of us originally used to show that the action 
for unparticles is identical to an action on anti de Sitter space[10]. This same model has been recently adopted[14] 
to obtain the hyperscaling violation exponent and anomalous dimension for the current. Within this approach, we 
obtain a power law in optical conductivity ~ uj~°‘ and the associated phase angle of The implication of this 

result for anomalous dimensions for the current[14] are discussed. 


II. UNPARTICLE 

Although one of us explicitly argued for the relevance of unparticles to the pseudo gap phase of the cuprates 
elsewhere[10], we summarize the key claim here for completeness. Numerous arguments have been presented for 
quantum critical points occurring in somewhere in the pseudo gap phase of the cuprates[4, 15]. However, precisely 
what are the propagating degrees of freedom has proven elusive to enumerate because the underlying system is strongly 
correlated. Ref. [10]’s arguments for the relevance of unparticles to the pseudo gap phase focused on the existence of 
Fermi arcs. Ref. [10] argued that because Fermi arcs require the presence of zeros of the single-particle Green function 
either at the terminus of the arc or at the backside of the arc, the excitations that give rise to them arise from a 
divergent self energy and hence are not connected adiabatically to a Fermi liqnid. In general the new excitations are 
some composites that are determined by collective phenomena. It is precisely composite excitations that conformal 
field theory (CFT) attempts to describe as the relevant degrees of freedom at critical points. Since any non-trivial 
infrared dynamics in strongly correlated electron systems must be controlled by a critical fixed point, scale invariance 
can be used to construct the form of the underlying propagator. It is this motivation that led Georgi[ll, 16] to forgo 
the full machinery of conformal field theory and just use scale invariance to construct the form of the relevant IR 
propagator. One of us put this proposal to use[10] in condensed matter systems because it is valid in any dimension 
and hence permits immediate progress beyond the restrictive ID Luttinger liquid paradigm advanced by Anderson[9]. 

Regardless of the explicit mathematical machinery (multi flavors, summation over masses, for example) used to 
construct unpaticles, unparticle stuff[ll, 16] is a scale invariant sector that emerges in the infrared (IR) of an effective 
field theory. Unparticles can be thought[17] to arise from a general field theory in which integrating out the ultra¬ 
violet (UV) sector generates new degrees of freedom at low energies, that is, UV-IR mixing. The signature of such 
mixing is an interpolating field that now resides in the low-energy sector [17] as in the charge 2e boson in the Hubbard 
model[10], the basic model used to describe the cuprates. This bosonic field mediates new electronic states at low 
energy that have no overlap with bare electrons and hence gives rise to zeros of the single-particle Green function, 
thereby signaling the breakdown of the Landau Fermi-liquid picture. We can understand the nature of the excitations 
that emerge by characterizing the new fixed point. In principle it is difficult to establish the existence of fixed points 
in strongly correlated systems. However, since all fixed points possess scale invariance, we can write the candidate 
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propagators down immediately. The propagator of a scalar unparticle in d + 1 Lorentzian space-time dimension can 
be written immediately from scale invariance as 

Gik) = —- (1) 

2sm(d;77r) - d.u 

with 

167r’^/^ T{du + k) 

[2TTYduY{du-l)T{2duy 

The unparticle scaling dimension djj must be greater than as required by the unitarity bound. When djj is not 
an integer, there are no poles in the propagator G{ko,k) when considering the complex variable ftp ■ Instead, there 
are singular points at ko = with a branch cut coming out from each point. In general, a single quasi-particle 
excitation is associated with a first order pole in the propagator. Hence, unparticle stuff is a theory with no single 
quasi-particle excitations. We introduce below two formalisms of unparticles that we use to calculate the optical 
conductivity. Indeed a well known example of unparticles in condensed matter systems is that of a Luttinger liquid. 
However, this is valid only in d = 1. A key result one of us also showed in the original paper on unparticles in 
condensed matter system[10] is that the continuous mass formulation of unparticles is equivalent to an action in an 
anti de Sitter (AdS)spacetime. That is, the summation over mass (or flavors) is equivalent to a summation over energy 
and hence the connection with the AdS radius. Hence, unparticles receive contributions from all energy scales. That 
is they represent the incoherence, the identifying feature of doped Mott insulators. Consequently, their relevance to 
cuprate physics is quite natural. 


A. Unparticle Effective Action 

Unless we have a specific conformal field theory model (e.g. see Ref [18, 19]), we do not know how the unparticle 
field couples to the 17(1) gauge field and how to write down the corresponding four-point correlation function. Both 
gauge coupling and four-point correlation functions are necessary to compute the optical conductivity. To this end, 
we write down a specific model with a Gaussian action, 

S = J (-27r)-^+i ^^(P)»G~^(p)^(p) (3) 

where (j){p) is a complex unparticle field in momentum space. This effective action can be used to correctly reproduce 
the two-point correlation function. One must keep in mind that this action is a particular model of unparticles 
in which the forms of higher-point correlation functions are specified by Wick’s theorem. An implication of this 
form of action is that an unparticle field is treated as if it is a free particle propagating with an algebraic form of 
momentum. We show how one can gauge this action and use it to obtain the conductivity in section HI. Nonetheless, 
the inherent particle content of this description leads to an underlying free form for the conductivity. This calculation 
is also instructive because it illustrates why the Luttinger liquid description of Anderson’s [9] cannot ultimately yield 
a power-law for the optical conductivity. 


B. Continuous Mass Formalism 

Unparticle stuff can be formulated as an integration over the mass of free fields [20, 21]. The general principle 
underlying this statement is that a scale invariant theory can be generated from a summation over the large number 
of flavors or from an integration over mass. Let (pi be a free scalar field of flavor i with mass mi, ranging from 0 to 
infinity, and let the density of flavors be The number of fields between to -I- is f{m'^)dm'^. An 

unparticle field can be defined as 

OO 

^uix) = Y^piix) = J dm"^ /(m^)())(x,m^). (4) 

* 0 

We will see that this definition is not the most general if unparticles serve to describe incoherent electronic states. As 
shown in Ref. [21], the unparticle two-point correlation function is given by 

OO 

G{p) = [ d'^^'^xe"'P'^{4>u{.x)(j)u{^)) = [ \ , ■ /^(^^)- (5) 

J J P — ^ 

0 
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For a choice of ~ 2 )/2^ the two-point correlation function is 


G(p) 


(—— le) 2 


^-d„ 


( 6 ) 


The interpretation of this formalism is that an unparticle held is an effective excitation (or effective held) of a contin¬ 
uum of free scalar helds with varying masses. The continuous mass formulation description was further generalized in 
Ref. [10] by introducing mass as another coordinate in a scalar held Lagrangian. The action can then be recast into 
Poincare coordinates of anti de Sitter space. Using the AdS/CFT correspondence, Ref. [10] calculated the two-point 
correlation function on the boundary theory and identihed the unparticle propagator. Ref. [14] showed that the 
hyperscaling violation exponent and anomalous dimension for the current can be generated by summing over the free 
energy in a multi-band system. Such an idea is completely equivalent to the continuous mass formulation and hence 
we will use the two approaches interchangeably here. Instead of calling it a multi-band system as in Ref. [14], we say 
that our system has a large number of havors of free helds. We apply the continuous mass formalism to show how one 
can get a power law in the conductivity in a system with a large number of havors (or multi-band system) in section 
IV. 


III. OPTICAL CONDUCTIVITY FROM UNPARTICLE EFFECTIVE ACTION 


In order to calculate a response function to an electromagnetic held, we need to know how the unparticle is coupled 
to a U(I) gauge held. To this end, we follow the approach for gauging a non-local Lagrangian with a Wilson line from 
Refs. [12, 22]. This approach was applied to scalar unparticles in Ref. [13] and to unfermion in Ref. [23]. To gauge 
the action (Eq. 3), one rewrites it in position space, 

S = J (l>'^{x)F{x-y)(f>(y) (7) 

where F{x — y) is a function resulting from converting the action to position space. The action can be gauged by 

V 

inserting the Wilson line, W{x,y) = exp{—ig f dw^A^), between {x) and (j>{y) with g being charge of unparticle. 

X 

The new action, 

S = J d'^^^xd‘^+^y (j)\x)F{x-y)W{x,y)(j>{y), ( 8 ) 

is now U{1) invariant. The vertex couplings of the unparticle and gauge helds can be obtained by taking derivatives 
of the action with respect to the helds. For the calculation of the optical conductivity, we need the vertex of two 
unparticle helds with one gauge held, 

S^S 

{Pid) SAt^{q)S(l)^p + q)S(l){p) 

= g{2p^^ + q^)F{p,q) (9) 

and the vertex of two unparticle helds with two gauge helds, 


g 2 F^'^(p, 91,92) = 


6^S 


SAf^{qi)6A‘'{q2)d(j>'^{p + qi + q2)S(j>{p) 

{2p + 92 )*^(2p -I- 292 + qi)^ ^ 
ql + 2{p + 92 ) • 9 i 
(2p-k9i)^(2p-k29i-^ 92 )“' 


=9^ 2p^^''F{p,qi +q2) + ' 2 ^' S"' ' i^iP^ li + 92) - d^(p,q2)) 


The function F{p, 9 ) is given by 




92 + 2(p + 9i) ■ 92 


^G-\p + q)-tG-Hp) 


{F{p,qi Fq 2 )-F{p, 91 )) . 


{p + gY — p^ 

These vertices satisfy the Ward-Takahashi identities, 

-Z 9 ^r''(p, 9 ) = G~^{p + q) - G~^{p) 


( 10 ) 


( 11 ) 


( 12 ) 
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and 


= r‘"(p + 91 , 92 ) - r‘"(p,92)- 


(13) 


The thermal version of the unparticle propagator and the vertex couplings, both of which are needed for the optical 
conductivity calculation, are shown in the Appendix B. The outline of the calculation is as follows. First, we compute 
a Fourier component of an imaginary time response function to an electromagnetic field from 




F2 SA^^ri{q)SA^^ri'{q') 


2 [A], 

A =0 


(14) 


where 2[A] = f D(j)D(f >^is a partition function of the action given in Eq. Bl. The subscripts, n and n', in 
the Fourier component of the response function above denote the dependence on the bosonic Matsubara frequency 
(wm = with m being an integer). In a system with translational invariance, the response function satisfies the 
relation 


Ar^''(Ti,Xi;T2,X2) = - T2,Xi - X2). 

Fourier transforming Eq. 15, one finds 

K:nM,q') = ^^"(9 + 9')W.oifr('?)- 


(15) 


(16) 


We use Eq. 16 to reduce from the Fourier component of K^'' {ti, xi]T 2 , X 2 ) to that of — T 2 ,Xi — X 2 ). After 

computing the response function, K^^{q), one obtains the conductivity in Matsubara space from 


q—^0 CJyj, 


(17) 


Here we take the the thermodynamic limit 9 —> 0. Finally, we perform an analytic continuation from the Matsubara 
to real frequency. Using Eq. 14, one can show that the response is of the form 

J j^^r((;^_„_„,(p,-9,-9')(9i'L-n-n'(P-9-9')</>m(p)) 

d^p d'^p' 


- E 


, y (27r)^ (27r)'^ 




m,m 

Expanding the four-point correlation function with Wick’s theorem and substituting 

\d 


{4L{p)(Pnip')) = - P')^m,nGjnip) 


(19) 


into Eq. 18, one finds 


K''{q)=g^TY, I ^C%.n(P,-'Z,9)G™(p) 

-S^T’E / ^J^^'L,niP^q)^m+n,-niP + q^-q)Gmip)Gm+n{P + q)- ( 20 ) 

The first and the second terms of Eq. 20 are known as the diagmagnetic and the paramagnetic terms, respectively. 
Summing over the Matubara frequency with the standard contour integration technique (see Appendix C), one finds 
the optical conductivity of the scalar unparticle (Eq. C31) is 

cr(zw„) = - du)(Joiiu}n) ( 21 ) 

where uo is the conductivity of a free massless scalar field (Eq. A 8 with m = 0). In d = 3, the prefactor is 2 — du- 
Liao[24] showed that the polarization is the same prefactor, 2 — d, times the polarization for the scalar particle 
held, thereby corroborating the lengthy calculation performed here. Note that, when djj < the real part of 

the conductivity becomes negative, which is not physical. Thus, the allowed value for du in a Gaussian model is 
< du < From the form of (Tq in Eqs. A9 and AlO, we see that cto has no dependence on du- Therefore, the 
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optical conductivity from the unparticle effective action does not have a fractional power law. In the context of the 
polarization for unparticles. 

At work here is the linear superposition that links the unparticle to the particle field, inherent in a Gaussian action 
for unparticles. Such a superposition implies that the statistics of the unparticle cannot deviate from that of the 
particles. Also at work here is the fact that vertex corrections within a Gaussian action kill the power law at the bare 
one-loop level. This can be seen as follows. The scaling dimension of the unparticle propagator is where 

[L] denotes a scaling in units of length. We know that the vertex couplings must satisfy the two Ward-Takahashi 
identities. So and have a scaling dimension of and respectively. As a result, the scaling 

dimension of the response function (Eq. 20) is and that of the conductivity (Eq. 17) is In other words, 

the anomalous length scale in the unparticle propagator, namely, anything that has a scale or its integer power, 

is cancelled by the vertices. This applies equally to the Luttinger liquid calculation [9]. 


IV. OPTICAL CONDUCTIVITY FROM THE CONTINUOUS MASS FORMALISM 

In this section, we apply the idea of the continuous mass formalism discussed in section IIB. We consider a model 
with large number of flavors of free fields or free particles. We sum over the conductivity of each flavor to obtain the 
total conductivity of this ensemble of free fields. In the continuous mass formulation of unparticles we explained in 
section IIB, there is no cutoff in the mass integration. Nonetheless as we will see, we take M to be a mass cutoff in 
the calculations below. 


A. Scalar Field 


We consider a model system with a large number of free complex scalar fields, 

S = ^ f dr f +mf\(l)i\'^), 


( 22 ) 


where the covariant derivative is defined in the usual way as — iciA^ and the summation is over all flavors. 

The mass and charge of the field of flavor i are given by rrii and e^, respectively. Let the density of flavors, p{m), 
to be a function of mass in the mass range from 0 to M. The number of fields in the mass range m to m dm is 
p{m)dm. Furthermore, we let the charge to be a function of mass, e(m). The forms of p{m) and e(m) are chosen in 
the same way as in Ref. [14]: 


p{m) 

m 

(23) 

e(m) 

O 

II 

(24) 


Here a and b are the power law scaling of the density of flavors and charges. As pointed out previously [14], a and b 
are related to the hyperscaling violation exponent and the anomalous dimension for the current, respectively. 

From Eq. A8 in appendix A, the optical conductivity of flavor i has the form 

(25) 

The total conductivity is a summation of conductivity of all flavors, 

= ^cr*(w) 


M 


= J dm p{m)e'^{m)f{uj,m,T). 


(26) 


We replace the summation by the integral over mass and the density function p(m). Rewriting the conductivity of a 
massive free scalar field (Eqs. A9 and AlO) in the form of Eq. 25 and scaling out the frequency, one obtains 


= ReCT*(a;) = ttpdS{uj) + Vi( —, -) 

UJ (jJ 

a^iuj) = ImCT*(a;) = — + e-a;'^"^/2( —, -). 

UJ UJ UJ 


(27) 


( 28 ) 
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The Drude weight p\) is given by 


Pd — 


2ejf3Sd 

d{2'KY 


' PI d+iN{eY{N^ 


1 ) 


(29) 


where ^{E) = E is the Bose-Eisnstein distribution, Sd = is a surface area of a unit (d-l)-sphere, and 

Ep = + mf. The functions /i and are given by 


/i 

h 


m T 

LV ’ UJ 

m T 

UJ ’ UJ 


TrSd 

(i2‘^(27r) 

Sd 

d{2TrY . 


-4 1- 


4m^ 


UJ^ 


dy- 




-P 


coth 


jduj 


(y2 + (^)2)i Vl_4(y2 + (^)2) 


coth 




V 


where 6{x) is the Heaviside function and P denotes the Cauchy principal integral. We consider a special case of d = 2 
and T = 0 (Eqs. A12 and A13). In this case, the Drude weight = 0 and the integral in /2 can be done analytically. 
One finds 


and 


/i 

/2 


1 

IGtt 


From Eq. 26, the real part of conductivity is 


i C \ 2 r / \ 

UJ 

i f \ 2 .f / \ 

=ej 2 ( —) 

UJ 


^1 -4 


1-4 - 


-4- + (l-4(-)2)ln 

UJ UJ 


1 - 2 ^ 


1 + 2 ^ 


M 


/ TTi 

dm p{m)e^{m)fi{ —) 

UJ 


(30) 

(31) 




i+2b 


^ ) 2“+2i>+3(a+26)(o+2&+2) M“+2i> 

_ 1 _ 1 t M\2 

4 


Poep ( 1 _ 1 ( 

4 l4(a+26) (a+2b+2)y u J ) 


if W < 2M 
if w > 2M 


(32) 


where a + 26 must be greater than 0 in order for the integral to converge. When uo/M < 2, cti has a fractional power 
law of However, when uj/M > 2, ui is approaching an asymptotic value with a power law of A plot of 

(Ti in the two-dimensional case at zero temperature is shown in Fig. 1(a). 

The imaginary part of conductivity can be evaluated approximately as 


M 


o' 2 (w) = luia(uj) = / dm p{m)e^(m)f 2 { —) 

J w 


(33) 


apoeo 


M 


5poeQ 


1677(0+26+1) u ) 2 “+ 2 ^+ 57 r(a+ 26 +l) V M 

Po^n M_ i‘f , , \ 9 A/f 

27r(a-|-26+l) uj n ^ > 2M, 


(^)“+2'> ifw<2M 


where a -I- 25 must be greater than 0 in order for the integral to converge. Unlike cti, when uj/M < 2, 172 depends on 
the terms ~ 0 ;“+^^ and ~ uj~^ . When uj/M > 2, CT 2 is proportional to uj~^. A plot of the approximated a 2 in the 
two-dimensional case at zero temperature is shown in Fig. 1(b). The real part of scalar conductivity in two dimension 
at finite temperature is 


cri(w) = 


Poep CP 


{ TrpD(j{uj) 

- {a+2b+2) (^f)<^Oth{^) ifw>2M, 


2o+26+3(p+2b)(a+26-|-2) M“+2i> 


coth(^) ifuj<2M 


(34) 




(a) Real part (b) Imaginary part 

FIG. 1. Real part and imaginary part of the total conductivity in 2D at zero temperature for the case a = 0.2 and b = 0.3 


where is the total Drude weight and is given by 

M 

PoellS 


Pd = 


27rM“+2'’ 


dm [ 


The imaginary part of the conductivity at finite temperature can be obtained from the Kramers-Kronig relationship. 
At finite temperature, cti is modified by a bosonic factor, coth(^). Consequently, for T ^ M, ai still has a power 
law of when oj < 2M. When T is about the same order of magnitude as M or larger, the power law in ui no 

longer exists. Hence, although we have obtained a fractional power law in the case of a scalar field, it is a power law 
with a positive power not negative power (a + 26 > 0) as required to explain the experiments. 


B. Drude Conductivity 


Motivated by the experiments [4], in which the real and imaginary parts have a negative power law, we apply the 
continuous mass formalism to the Drude form of the conductivity. The system consists of large number of flavors of 
free non-relativistic fermions in a disordered potential with no interaction among different flavors of fermions. The 
conductivity of one flavor is simply a Drude conductivity [25], 


a\u}) = 


UiefTi 


1 


(35) 


mi 1 — iuTi 

where mi, e^, rii, and Ti are mass, charge, density, and relaxation time of fermion of flavor i. As in the case of a 
scalar field, we take the density of flavors, p{m), charge e(m), and relaxation time, rim), all to be functions of mass. 
The mass ranges from 0 to M. The charge density in the mass range of to to m + dm is dn(m) = p(m)dm. So the 
contributions to the conductivity from the flavors in the range to to to + dm is 

p{m)e‘^(m)T{m) 1 


da = 


TO 


1 — iujT{m) 


dm. 


The conductivity of all flavors is 


M 


a{uj) = J 


p{m)e‘^{m)T{m) 


1 


TO 


1 — i(jjT(m) 


dm. 


(36) 


(37) 


We choose the form of the charge, density, and relaxation time in the same way in the scalar field case, namely. 


p{m) = po- 
e{m) = eo 




T{m) = ro^. 


(38) 

(39) 

(40) 
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The conductivity becomes 


a{uj) = 


PoepTo 

]\^a+2b-\-c 


M 


dm 


m 


(i-\-2b-\-c —2 


1 icOTQ 


Changing the integration variable to a; = wrp ^, one obtains 


r(w) = 


PoCp 


1 


U)To 


cM 


w(a;To)‘ 


dx 


1 — ix 


(41) 


(42) 


Define a = 1 . When a < 0, the integral does not converge. This means it is not possible to obtain a positive 

power law in the optical conductivity from this model. When 0 < a < 1, the integral has a simple form in the high 
freqeuency limit WTp 1: 


cr(w) 


PoSq 


Tre 2 


-U) 


sinTra 


(43) 


OO 

Since ReCT(a;) falls off slower than the total spectral weight, f Reo'(uj)du! is infinite. So we need to define a high 

0 

frequency cutoff, ujc, to have a finite spectral weight.^ Here we see that by choosing the value of a, b, and c such that 
0 < a < 1, the conductivity can have a fractional power law with a power between -1 and 0 and a phase of 

As a concrete example, let us try to reproduce the power law of — | seen in the optimally doped cuprates above 
the superconducting dome [4]. We choose = — | with c = 1, so a = |. The conductivity is 


cr(w) = 


PoepTj* 


M 


dx 


(jJ3 


1 — ix 


PoepTj* 


M 


ln(l + i(Toa;)3) + e 3 " ln(l — e" {tooj)3) — e" ln(l — e " (Toa;)^) 


(44) 


where the branch cut of the logarithm is chosen to be along the negative real axis. The cutoff is set at some high 
frequency loc- Taking the limit of tquj —)■ 00 , we obtain 

(t(w) = ^(v^ + 3i)7r ^°^°^° . (45) 

3 Muji 

We see that the power-law exponent is — | with a phase of 60° as is seen experimentally. While there are ways 
of generating fractional exponents[26, 27] in the optical conductivity, none yield the exponent of —2/3. The recent 
claim that gravitational crystals[28] produce the conductivity of the cuprates was not borne out by subsequent 

calculations [1-3]. Consequently, the results obtained here are the first to yield the —2/3 exponent from a strong 
coupling formulation. Our result is general, following strictly from the continuous mass formalism and the Drude 
form for the conductivity. Hence, what is needed of a microscopic theory is one in which the optical conductivity at 
each frequency receives contributions from all energy scales. That is, incoherence is at the heart of the power law 
optical conductivity in strongly correlated electron systems. 

When a > 1, the conductivity still has a power law of in the real part and there is no need for the high 
frequency cutoff. For example, we choose | with c = 1, so a = |. The conductivity in the limit of high 

frequency is 


cr(w) 


Poeo,5i 

M^oj 


1.65163- 5.0832t^ 

Tq to 5 


(46) 


The real part of the conductivity has a power of — whereas the imaginary part has an extra term ~ cc 


^ Even though the spectral weight of one flavor is finite, 



1 


l + cjVP 


duj 


TTUief 

2mi 


the integral of the total spectral weight, f diverges when 0 < a < 1 and c is a positive number. Here c must be 

0 

positive, otherwise the real part of conductivity would be negative. We still need to define a cutoff, order for the total spectral 

weight calculated this way to be finite. 
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V. DISCUSSION AND IMPLICATIONS 

The result that there is no fractional power law from a Gaussian form of unparticle action is not a surprise as 
follows from the scaling form of the vertex. On the other hand, other response functions that do not have this 
dimensional cancellation from vertices should exhibit a non-integer power that depends on du. One example of this 
kind of response function is the density-density correlation {(j)’’' {x)(f>{x)(j)^{y)(f>{y)). 

From this reasoning, one might think it is possible to obtain a power law if we know a way to embed an anomalous 
dimension for the current into the effective action. The present approach we used to gauge the action using Wilson 
lines does not allow one to do so. If we can find another way to write down a gauged effective action with the 
anomalous dimension for the current included, one should be able to get a power law that depends on the anomalous 
dimension. However, from the result of the mass summation over the two-dimensional scalar field at T = 0, we 
find that the exponent is a -I- 2b, where a is related to a hyperscaling violation exponent (or djj) and b is related to 
anomalous dimension for the current. Clearly, even when & = 0, it is still possible to have a power law. This signifies 
that one can nevertheless obtain a power law in the optical conductivity without the inclusion of the anomalous 
dimension for the current. 

The fact that the conductivity of the Gaussian unparticle is the same as that of massless scalar particle up to a 
constant multiplication factor is interesting. The technical reason of this result is due to an exact cancellation of 
discontinuities across branch cuts of diamagnetic and paramagnetic terms. The response functions that do not have 
this kind of cancellation, such as density-density correlations, should acquire non-trivial dependence on the anomalous 
dimension, djj, unlike the case of the optical conductivity. 

We want to mention several points about the scalar field summation results. First, in d = 2 and T = 0 case, ai has 
a positive power law. This is due to the fact that we sum over a collisionless contribution to the conductivity (see [29] 
for an explanation of each part of the conductivity). Should we add a momentum relaxation process, for example, 
dissipation and interactions, and then sum over this contribution to the conductivity we would get a negative power 
law. Phenomenologically, this part of the conductivity should have a Drude form. We perform a similar calculation 
on the Drude conductivity in the section IV B. We find that, indeed, the conductivity in that case has a negative 
fractional power law. The second point is that ci in the case of d = 2 and T = 0 has a fractional power law, whereas 
(72, in addition to a term ~ has a term that goes like ~ cv~^. At low frequency, CT 2 is dominated by U!~^ 

behavior. At higher frequency, the power law of could appear in the frequency range M < oj < 2M, if the 

exponent a -I- 26 is large enough. This result suggests that in a system in which cti has a positive fractional power law, 
CT 2 may change from having an integer power law at low frequency to a fractional power law at higher frequency. In 
such cases, the power law disappears at finite temperature due to a bosonic temperature factor, coth(l^). As long as 
the temperature is kept to be much lower than M, there is still a power law in cti. When the temperature is raised 
to be about the same order of magnitude as M or higher, the power law disappears. 

The key result of this paper is that continuous mass formalism coupled with the Drude conductivity yields a power 
law with exponents between -I and 0 which are in agreement with the exponents seen in the cuprates. The key point 
of the continuous mass formalism is that summation over mass is summation over energy. Hence, the power law 
found here arises from the incoherent part of the spectrum, that is the part that has no particle content, hence the 
description in terms of unparticles. As pointed out in section IV B, one needs to set a high frequency cutoff in order to 
have a hnite spectral weight. While the continuous mass formalism works well as a phenomenological tool to explain 
an appearance of power laws, it cannot produce other features that require microscopic details. Our results suggest 
that a microscocpic theory that can explain fractional power laws should be decomposable into a large number of free 
fields or free particles with varying masses and possibly varying charges. Still, one should keep in mind that there 
could be other models that predict a power law. We should also note the result when the exponents are less than 
-1. (Ti has a fractional power law oj~°‘ at high frequency, but <72 behaves like Au!~°‘ + Blo~^ where A and B are real 
constant. 

To conclude, we studied the optical conductivity in two systems: a system of unparticles with a Gaussian form for 
the action and a system with a large number of flavors of free fields. We found that the conductivity of the Gaussian 
unparticle does not have a fractional power law that depends on du. However, using the continuous mass formalism 
on the system with a large number of free helds, we obtained a positive fractional power law in the case of free scalar 
field and negative fractional power law in the case of Drude conductivity. That the combination of exponents here 
can be constrained by other transport measurements (Hall angle and DG conductivity) on the cuprates will be the 
subject of a forthcoming publication. The latter result points to two possibilities for a microscopic model: either an 
anomalous dimension is present which indicates some kind of multi-band system or the anomalous dimension arises 
from UV-IR mixing. This has been pointed out extensively in the Hubbard model[10]. Hence, perhaps it might be 
able to construct the power-law conductivity from a theory in which the upper Hubbard is integrated out. Ref. [30] 
relied on the non-crossing approximation and hence any effort forward must be based on new techniques that can 
incorporate the vertex corrections when a particle-picture breaks down. 
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Appendix A: Optical Conductivity of a Free Massive Scalar Field 


The optical conductivity of the free massive scalar field has been studied by many authors (for exmaple, see Refs. 
[29, 31]). Since the scalar field conductivity is used in section IV and appendix G, we briefly review the basic 
formulation and the main results in this appendix. The optical conductivity can be calculated from Eq. 20. The 
propagator of the massive scalar field in Matsubara space is 


G„(fc) = 


1 


+ et 


(Al) 


where the bosonic Matsubara frequency (n being an integer), k is a d-dimensional momentum, Sk = 

The vertices in this ( 


+ rn?, and m is the mass of the scalar field. The vertices in this case are 




(A2) 

(A3) 


For the expression is a Matsubara frequency and p*, for t = I to d, is an ith component of d-dimensional 

momentum. Substituting Eqs. Al, A2, and A3 into Eq. 20, one finds the longitudinal conductivity (the ii component 
with z 7 ^ 0) in the limit g —)■ 0 to be 


V 

UJr, 


a{iuJn) = —T 


(Yp 

(27r)^ 


1 


2 pf 




Inserting I = ^ into the diamagnetic term and integrating by parts, one finds 

f Yp 1 I 


cr(za;„) = - T 

The boundary term ^ f 


2 

(27r)^^* V « + 


(A4) 


(A5) 


P 


^2 +52 |^-_oo vanishes in the dimensional or lattice regularization [29]. We can 

also see this clearly in the case of a lattice system with periodic boundary conditions. In this case, pi is replaced by 
sin Pi and ±00 becomes the momenta at the edge of the Brillouin zone = ±7r. As a result, the boundary term vanishes 
(see appendix B of Ref. [32]). Had we not done the integration by parts, the optical conductivity would have a 
zero-frequency delta function peak at zero temperature. We implicitly set the chemical potential of the scalar field to 
0. According to the Bose-Eisenstein distribtution, there are no particles or antiparticles present at zero temperature. 
This means that the dc conductivity should vanish and thus the zero-frequency delta function is not physical. At 
finite temperature, the mixture of particles and antiparticles is excited into the system and so the zero-frequency delta 
function should appear. 

We perform Matsubara a summation explicitly using contour integration. Let us first consider the summation on 
the diamagnetic term. Converting the summation into a contour integral with the contour shown in Fig. 2(a), we 
have 


T 


1 




\^m ' 


r dz g{z) 
Ic2m ( 02 _^ 2 ) 2 - 


(A6) 


with g{x) = I coth(^). Upon deforming the contour C to infinity, the contour integral turns into a summation of 


the residues at 


P’ 




'(jj‘^ _L £-2\2 
' ^pJ 


E 

Zi=±e„ 


ReSz 


5(2 


{z^-elY 


I 


— + 1 ) + i^dYp)^ 
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(a) Diamagnetic term (b) Paramagnetic term 

FIG. 2. The contours used in the Matsnbara snmmation of the diamagnetic and paramagnetic terms. 


where N{x) = is the Bose-Einstein distribution. The similar computation, with the contour shown in Fig. 2(b), 

can be done on the paramagnetic term to show that 




1 

+ep)((Wm 




gjsp) f 1 _ 1 \ 

iujfiSp ^icOji ‘2sp T ‘^^p) 


(A7) 


Combining the results for the paramagnetic and the diamagnetic terms, one finds the conductivity in Matsubara space 
to be 


/ \ r. 2 /" 2 

aiiuJn) = 2 g J 


l^Niep)iNiep) + l)^'-^ 


^nP- 


1 




Pefroming analytic continuation iwn —)■ w + tO’*', the conductivity turns out to be 

g^T^Sd 


ai{uj) = Re(T(w) = irpoSiuj) + 

. N j /X Pd , g'^Sd 


d2d{2TrY 


0 (a; 2 -( 2 m) 2 )a;"- 2 (l- 

r, 




where 9{x) is the Heaviside function, P denotes the Cauchy principal integral, Sd = 
(d-l)-sphere and the Drude weight po is given by 


1 

2 £p — iio 

!)<i/2coth(^) (A9) 

(AlO) 

d 

the surface area of a unit 

^ \ 2 ) 



OO 

V/3Si j+iiV(ep)(iV(e,) + l) 
- 4 - 

0 


(All) 


Note that at T = 0, p/j = 0. So the zero-frequency delta function vanishes at T = 0 not at finite temperature. There 
is a gap of 2m in the real part of conductivity as can be seen from the Heaviside function in Eq. A9 (also see Fig. 
3(a)). This gap is the minimum energy required to create a particle-antiparticle pair in the system. In the special 
case of d = 2 and T = 0, the optical conductivity can be calculated analytically. 


CTi(w) 

CT 2 (w) 






4to 


(2m)" 


)ln 


oj — 2m 
(jj + 2m 


)■ 


(A12) 

(A13) 


The plots of the real and imaginary parts of conductivity in this case are shown in Fig. 3. 
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(a) Real part (b) Imaginary part 

FIG. 3. Conductivity of massive scalar field in d = 2 at zero temperature. 


Appendix B: Unparticle in Matsubara Space 


The U(l) invariant unparticle effective action in d+1 Matsubara space is given by 

+ / ^y^^ln+niP + <l)(t>m{p)Af,,n{q)g^';n,n{P,Q) 

+ ^ ^^ '^m+ni+«2(P+gl +g2)<?^m(p)^p,ni(gi)^^.,n2(g2)yr^'^ni,n2(P^gl^g2)> (^1) 

m,ni,n 2 ' 


where subscripts of the helds, the propagators, and the vertices denote the dependence on the bosonic Matsubara 
frequency (a;„ = with n being an integer). The propagator is 


G„(fc) 


^du _ ^ _ 

2sin(d,77r) +ujl)^-du ' 


The two unparticles one gauge field vertex is 

q) = ( 2 p'" + q) 

and the two unparticles two gauge field vertex is 


(B2) 


(B3) 


^ (2p + g2)"(2p + 2g2 + gi)^ 

9l + 2(p + ( 72 ) • qi + + 2{uJrn + 

_l_ _ (2p + gi)^(2p + 2gi + ( 72 )'^ _ 

+ 2 (p + qi) ■ q 2 + + 2 (uj m H“ ^n\ )^n 2 


m,ni+n2 {.Pi Ql 
(‘^m,ni+n 2 {Pi Ql 


Q 2 ') ^ m ,712 ( p > 92 )) 

“(“^72) ‘ 7 'm,ni (P; 


(B4) 


where the function J^{p,q) is 


^m,n{Pi Q) 


Gj+niP + g) - Gj{p) 

q'^ + 2p - q + ujl + 2uJmUJ„ 


(B5) 


From the expression is Matsubara frequency and p*, for f = 1 to d, is an ith component of d-dimensional 

momentum. 
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Appendix C: Calculation of Optical Conductivity from Unparticle Effective Action 


From Eq. 20, we refer to the first term as the diamanetic term, „(<?)) s-nd the second term as the paramagnetic 
term, Substituting Eqs. B2, B3, and B4 into these two terms, one obtains 






{2p + qY{2p + q)f^ - du) 

q"^ + 2p- q +col + 

{2p + qY{2p + qY 

(g 2 + 2 p • g + w 2 + 2 wmW „)2 


( 2 p + g)‘^( 2 p + g)^ / {p + g)^ + {uo^n + ^«) 


(g2 + 2p • g + w2 + 2 uJmO}nf 


d + 1 j N 

2 \ — - du 


p 


and 


^ ( 2 p^ + g^)( 2 p- + g-) 

para,nW 9 2^ J ^2^)^ \ (^q2 + 2p ■ q + Col + 2cO^COnY 


1 {2pt^ + qYi2p'^ + qn f ip + qY + {cOm+i^ny 

2{q‘^ + 2p-q + ujl, + 2LOmUJnY \ 


1 i2p^ + qYi^p''+ qY 


P^+ujI, 


2 I 2 

p +‘^t 


— du 


— du \ 


2 {q‘^ 2p ' quj^ 2u;^a;^)2 y (p + q^ + (cjm + u^ny 
In obtaining the result for the diamagnetic term, we interpret J^m,o(p,0) 


1 • ~r~ t \ ^ / 2 \ 


P.\^-dv 


q—^0,n—^0 


2^2—p2_^^2 


yd 1 ,\/9 9 \ 

= - du)ip +W„) 2 


We split K^ia into four terms and K^ara into three terms: 

T^p.u , _ _9 2 d+1 _ \ f d'^p {2p + qY{2p + qY _ 1 

d,a2,nW- 2 (27r)<ig2 + 2p-g + a;2+2a;„cc„<+p2 

K^'' («) = -2g2TV f _ {2p + qYi2p + qY _ 

d*a3,nW 9 2^j ( 2 ^)d(g 2 + 2 p.g + ca 2 + 2 w™a ;„)2 

(q) = 2q^TY^ [ _ ( 2 p + g)‘^( 2 p + g)^ _ 

d^a4,nW 9 2^j i2Y'^(^q^^2p-q + col + 2oo^ooY 


d‘^p {2p + qY{2p + q)f^ ( {p + qY + Ym + 0On)‘ 

p'^+ujY 


- —<i[/ 


K^aralAd) = -KdL3,n{9) 

1 


K^ara2,ni9) = - l^KdiaA^nid) 

n- ^tV f (2p + g)^(2p + g)^ / 

paras, nyq) 9 (^2nY + 2p ■ q + CO^ + 2cOmCoY^ \{p + qY + Ym + 


i-dr. 


(Cl) 


(C2) 


(C3) 

(C4) 

(C5) 

(C 6 ) 

(C7) 

(C 8 ) 

(C9) 

(CIO) 


We show here the Matsubara summation using the contour integral calculation technique on these seven terms. 
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(a) Kdial (b) ^dia2 (c) t^diaS 

FIG. 4. The contours used in the summation of three diamagnetic terms. 


1 « ^dial 


The ii component of Kdiai in the —)■ 0 limit is 




2 &’’)-<+p‘' 

Inserting a factor 1 = ^ into the integrand and then performing an integration by parts, we have 

'd+1 , f d'^P ‘^P'i 


Klaijq ^ 0 ) = 2 g^ 


-du)Tj2 


(27r)<^ +p 2)2 


Carrying out the Matsubara summation using the the contour shown in Fig. 4(a), one obtains 

1 1 






= ^E 




= ^^(p)(^(p) + 1 ) + ^9{p), 


where g{x) = ^ coth(^) and N(x) = is the Bose-Einstein distribution. The result is 

'd+1 \ r 


Klai.niq ^ 0 ) = 


<^u) + + (Cll) 


2 . Kdia2 


We rewrite the ii components of Eq. C5 as 


K, 


dia2, 


M = 


-w 

UJn 


--<b)TY. 


d'^P i ,x2 1 1 

(27r)'* ZWm +i(d + ^) 


where 5 is defined as 

2w„(5 = q^ + 2p-q. (C12) 

We will not take the limit 5 —)■ 0 (or q —>■ 0) right away, but we will do so after combining the contributions from Eq. 
C5 - C7. The reason for this is given below. We convert the summation into the contour integral with its contour 
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shown in Fig. 4(b): 

+ z{S +^)u;l+p^-~ Jc + z(^f + J) z 2 _p 2 • 

In addition to the poles of g{z), there are poles at ±p and —*(^ + S). Had we set <5 = 0, there would be a pole at 
—This is problematic. The reason is that when n is odd, the integrand has a simple pole at —i^- But when 
n is even, the integrand has a double pole at because ^ would correspond to another Matsubara frequency 

uji with I = ^. This problem is related to the fact that the summation T ^ |,iJn - ^2 ^^2 diverges when n is even. 

Hence, we need to keep S finite and then take the limit <5 —)■ 0 after we combine Eqs. C5 - C7. Deforming the contour 
to infinty and then summing the residues at ±p and —+15), one finds 


K. 


dia2. 


^(9) 



Wn 



dfip 

(27r)‘^ 


{2f + q^f 


i COt(f ((5+ f)) _ g(p) _ g{-p) \ 

2p^ + {5+ ^ 2p{p + i{5 + ^)) ^ 2p{p - i{6 +^))J' 

(C13) 


3. KdiaS 


We rewrite Eq. C 6 using the definition of S given in Eq. C12 and then perform the Matsubara summation with 
the contour given in Fig. 4(c): 


K. 


3.n(9) 


d‘^p 


dia‘i,n\'il o, ,2 ^ 2 !^ / (27r)<^ 

f d’^P 


i\2 


{2Y + qY 


1 


{iuJm + + ‘^)Y 


2 I P / r- 


8 ujI J { 27 rY 


i2Y + YrcsY YiS+^) . 


(C14) 


4. KdiaA 

We rewrite Eq. C7 with the definition of 6 given in Eq. C12 and then convert it into the contour integral with the 
contour shown in Eig. 5: 


KUn(q) = - 


g^ f d’^p 

2ujI J {2TrY 

g^ f d<^p 


J\2rp \ ' 


i2Y+YrT^ 


1 


1 


{iuJm + *(<5 + ^)Y 


= -^ + dz 


2colJ i2nY 


2TTiJc iz + i{6+'^)Y 


p^ + 2uJn5 - {iiOm + i^nY 
p2 + 2 w„(5 - (z + iiOn) 

p 2 — 


d +1 J 

2 \ -^-du 


(C15) 


In addition to the pole at —i{6 + j), there are branch points at ±p and _|_ 2uj„S — iu>n- Thus, we need to 

consistently choose the Riemann surface that will be used in the contour integration. We split up the term with 
power — du in the integrand as follows, 


p2 + 2uJnd - {z + iijJnY 
p 2 — z^ 


--dr 


{z -b icOnY - {P^ + 2.UJn6) 


z 2 — p 2 


{Z + ioJn - + 2uJn6) 2 _|_ 2uJnd) 2 


{z-pY 


’{z+pY- 


The constraint on the choice of the Riemann surface is the condition that the contour integral is real. In other words, 

dTi ^ sn ^ 1 to De 


we require the residue at z* = itOm of Eq. CIS, that is 
real. 




We choose branch cuts as in Eig. 5. One choice of the definitions of the phase angles that are satisHed by the reality 
constraint mentioned above is listed below. 
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FIG. 5. The contour of Kdiai in the case of odd n and small 5. The paths that go along the branch cuts are Ai to A 4 , and the 
paths that go around the branch points are Ci to C4,. 


For the term (z + , the dehnition of the phase angle is —tt < 0i < tt. 

For the term {z — , the dehnition of the phase angle is 0 < 6*2 < 27r. 

For the term {z + the dehnition of the phase angle is —tt < 0^ < tt. 

Lastly, the phase angle is 0 < 6*4 < 27r in the term (z — . 

The contour integral can be split into the contribution from the residue at —i{6 + y), the contribution from the 
discontinuity across the branch cuts {AiS in Fig. 5), and the contribution from the small circles around the branch 
points (Ci’s in Fig. 5). The residue contribution to the Matsubara summation is 


— Res[ 


(z + ziS+^W 


p^ + 2w„(5 - iz + iUTn)" 

p^ — 


- dr) 




^ ^cot{l{6+^)) 
(C16) 


Hence, the residue contribution to KdiaA is 


^^diaA^res^n 


( 9 ) 


9 ^ 

2 w 2 


d’^p 

{ 2 ttY 


(2/ + Yf 





du) 


cot(f(5+^))\ 
p2 + (5+^)2/ 


(C17) 


For the contributions from the discontinuity across the brach cuts and the small circles around the branch points, we 
take the limit q —> 0 or 5 —)■ 0 to simplify the calculation. Let us consider the contribution from the discontinuities 
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across the branch cuts (the integral along Ai to A 4 ): 


—p—s 

= 7 ^ [ dz 

2tti J 




gjz) iz + iLOr.-p)^-‘^^{z + ico„+p)^-^^ 

{z + i‘^y (_2 — 

g{z) {z + iu}n-p)^~‘^"{z + iu}n+p)^~‘^" 


{z+ + p) “*2 ^ (z_ + p) “*2 ^ 


27ri J (z + 1^)2 

p+e 


(z+p)'^ 2 ^ 


— l<jJn —p — e 


-—di) 


+ 


2'Ki 


dz 


— HjJji — OO 
— icjji+oo 


+ — / dz 


27ri 




gjz) _ (z + ia;„ - p) 2 _ 

(z + i^)2 +p)^~‘^'^ 

gjz) (z + ^u:r,+p)^^-du 

(z + i ^)2 — p'^^-du +p)^~‘^'^ 


(z_|_ — p) ”^2 ^ (Z_ — p) "^2 ^ 


/ • r/rr / • X d+'i 

(z+ + IW„ + p) 2 - [z_ + lOJn + p) 2 


(Z+ + - p) 2 - (2;_ + iwn - p)' 


— dr 


■t - d,u 


where z± = z ± try with ry —> O'*". Using the definitions of the angles 6 i defined above, we find that the discontinuities 
across the four cuts are 


1 1 

(z+ + p)^-'^’^ ~ (z_ + p)^-'^’^ 

1 1 

(z+ + itUn + p) ^ djj _ _|_ _|_ 2 ‘^u 

(z+ + iw„ -p)^- (z_ + iujn -p)^ 


—2zsin7r(d^ — du) 

2isin7r(d^ - 

|z — p\'^~’^'J 

2 i sin7r( ^^ ^ - d( 7 )|z + zw„ 

— 2 zsin 7 r(^^^i-|z + iw„ — p\^~'^'^. 


(C18) 

(C19) 

(C20) 

(C21) 


Substiuting the four discontinuities into the integral along Ai to A 4 and then changing variables so that the limit of 
integration is from £ to 00 yields 


simr(d^ - du) f,f 9{z + p) (-z + za;„ - 2p)(-z + zw„) “* 2 ^ 

^ J [(-z + z^-p)^ (^_z-2p)^-<^-z^-<^- 

gjz+p) jz + jz + iuJn + 2p)^-^^ 

(z + z^+yz)2 z^~‘^’^ {z + 2p)^~‘^’^ 

g{z+p) i-z-2p)^-'^^z^-‘‘^ 

(-Z - z^ -p)2 _ 2p)^-du(-z - 

_ g{z+p) _ z;^-rfv(z + 2p)^-^^ ^ 

(z - z^ +yz)2 'y 


On the third and fourth lines, we use the fact that g(z + iuJn) = \ coth(^^^^ 2 h!f£zil^ gj.g^ ^j^g 

third lines, we use g(—z) = —g{z) . To further simplify the integral, we pull factors of —1 out from the terms with 
power of — du and -1 in front of z to make the coefficients in front of z all be +1. At the same time, we redefine 
all the angles 9i to be evaluated by the power — du to be in one range, i.e. —tt < 9i < tt. These operations 
amount to multiplying each line by a factor d^'^j'^i^-du) (subscript j means line j). We find that mi = 1, m 2 = 0, 
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m 3 = 1, and 1714 = —2. The result is 


simrC^-du) f ( g{z + p) {z - iujn + 2p)''t^ {z - 

^ J \{z-i‘t+pr {z + 2p)^-<^-z^-<^- 

gjz+p) {Z + jz + iuJn +2p)^-‘‘" 

{z + i‘^-\-pY z^~'^’^ {z+ 2p)^~‘^’z 

_ g{z+p) _ {z + 2p)^-‘^^z^-‘^^ 

{z + + py [z + iujn + 2p)^~‘^’^ {z + 

_ g{z+p) _ z^-‘^^{z + 2p)^-‘^^ 

{z — +py (z — iuJn)^~^'^ (z — iuJn + 2p)^~‘^’^ 


Note that the first and the third terms are complex conjugates of the second and the fourth terms, respectively, and 
hence the result of the Matsubara summation is real as we suspected. Finally, let us consider the contribution from 
the integrals along the small circles around the branch points (Ci to C 4 ): 


1 

_ j 

f 9{z) 1 

((z + zw„)^ -p'^\ 

2TTi J 

c. {z + i'^y^ 

[ z^^-p^ ) 


(C22) 


For the curves Ci and C 2 , we let z = =FP + where — refers to Ci and + refers to € 2 - The integral is along the 
clockwise direction; hence it needs to be multiplied by a factor of —1. The result is 


£ 2 


- -\-dif 


Sf 


27r 


d6e 


id 


g{zs^p + £e*®) ( (=f2p + iuJn + £e*®)(fa;„ + £e*®) 


(=FP + f^ +£e*®)2 \ 


(t 2 p- 


■ £e*^)e*^ 


— di) 


(C23) 


In the case of curve Cz and C 4 , we let z = =FP — + £e*^, where — is for Cz and + is for C 4 . The result is 


£ 2 


- —djj 


27r 


dOe 


(t 2 p + £e*®)e*® 


10 g{TP + ge^ ) 

(=FP — + £e*®)2 (t2p — iw„ + ee^^){—iujn + £e*®) 


— du 


(C24) 


Eqs. C23 and C24 vanish in the limit £ —> 0 when the exponents of g^+du ^j-g greater than zero, that 

is, when < du < From the unitarity bound, du > If we now also require that du < we hnd that 
there will be no contributions from the Cj’s. 

Hence, the contribution from the branch cuts (both H^’s and Cds) to KdiaA in the limit g —?■ 0 is 


K, 


diaA.cut.n 


{q^0) = -g^ 


vSmTTl 


i^-du) 


2tt uj"i 


d’^p 

(27r)' 


KV 


.^2 


dz 


( g{z+p) {z-iU}n+2pYy {z - 

\{z-i^+py (z + 2p)^-^^z^-''^ 

g{z+p) ( 2 : + fa^n)'^~^^( 2 : + fa^n + 2 p)^-‘^^ 

{z + i‘^+pY {z+ 2p)^~‘^^ 

_ g{z+p) _ (z + 2 p)^-'^^z^-^^ 

(z + +py (z + iuj^ + 2p)‘^~‘^’z {z + 

_ g{z+p) _ z^-‘^"{z + 2p)^-‘^^ 

(z — + py (z — ia;„)^“'^"(z — za;„ + 2 p)^“‘^^ 


(C25) 



20 


5. 




Using Eqs. C 8 , C9, C14, C25, and C17, we find 

Kparal,ni<l) = ^ J (V + 9*)^ CSC^ 

^ + 9T ( - f csc^ {liS + f)) - 


(C26) 


du) 


cot{l{S+‘^))\ 

p^ + is+‘^rj 


(C27) 


and 


d^para2^cut^ni.^ ^ 9 


2 sinTT- du) f d‘^p 


dTTW^ 


(27r)< 


{2f + / dz 


(z-i^ +p )2 


g(z+p) (z - + 2 p)(z - iWn) ■'J^ 

, d+1 


(z + 2 p) ”^2 ^ '^c/z:”^ 2 ^ 


du) 


gjz+p) jz + iujr,)''^" jz + iuj„ + 2p)''t" ip. 


(z + i‘^ +p)‘^ 
g{z + p) 


z”^2^ (z + 2p)'^2^ 
(z + 


(z + +p )2 ^2; + IWn + 2p) “*2^ '^'^(z + fw„) “*2^ 

g{z+p) z^~‘^'^ {z + 2p)^~'^^ 

(Z - +p)2 + 2p)^-‘i^ 






(C28) 


For Kparas, we rewrite Eq. CIO in terms of 5 defined in Eq. C12, 

KparasAd) =^j ^ ^ + i (d + )) 


2 I 2 


— d[/ 


^))^ + 2w„(5 - (fwm+*w„)2y 

The calculation proceeds in the same manner as in the case of K^iai- We simply state the results here. The contribution 
from the residue is 


-a t \ 9^ f d‘^p j U2f /? 2 /^/d/c ,d+l cot(|((f+^)) 

para3,res,nid) - ^ J J^i^P + 9 ) “ 4 ^SC 2 + T 7 ^^ ^ p2 + (J + <^)2 


(C29) 


and the contribution form the branch cut (with < du < is 


7s:: 


para3,cut,n 


(q 0 ) =g‘ 


,sin7r(^ - du) 
dTrw? 


d^^p 

(27' 


K27 + g 0 " / 


/ 9 ( 2 :+p) (z - ia;„ + 2p) (z - ia;„) 

l^(z-f57+p)2 

9{z+p) jz + iuJn)'^~‘^'^ jz + iuJn + 2p)^-‘‘^ 

{z + i‘^+p)‘^ zT^“'^"(z + 

_ 9 {z+p) _ (z + 2p)7^-‘^^z7^-^^ 

(z + +p)^ (z + ia;„ + 2 p)T^“'^'^(z + 

_ 9 ( 2+7 _ 2:7^-^c/(z + 2p)7^-^^ 

(z — + p )2 ( 2 ; — (z — iujn + 2p)^~^^ 


(C30) 
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6. Total Response Function and Conductivity 

The contributions from the residues at —+d) vanish for both diamagnetic and paramagnetic terms, separately. 
To check this, we add the contributions from the residue at + S) from Kdia 2 , Kdias, and Kdiai,res, 


1 ,d+l COt(§(d+t)) /? 


(27r) 




2 su,d“' 




OJn 


+ ^esc^(^l(d+^))+—( 


1 ,d+l 


20Jn 


- du) 


cot(§((5 + ^)) 
p2 + (,5+5f)2 


= 0 , 


and add the contrubutions from the residue at —+ 5) from Kparaii Kpara 2 i and Kparaz-, 


d‘^p 

{2TrY 


{2f + q^) 






Wn 


16a; 2 


CSC -(<5 + —) 


1 ^d+1 j ^ COt(§((5 + ^)) 


16a; 2 


Ur, 




1 ^d+l 


- du) 


cot(f(d+^t)) 

p 2 + (5+!f)2 


= 0 . 


We can safely take the limit g —?• 0 and d —>■ 0. Futhermore, the contributions from the branch cuts, when adding 
both paramagnetic {Kpara 2 and Kparas) and diamagnetic terms {KdiaA), vanish, i.e. Eq. C25 + Eq. C28 + Eq. C30 
is equal to zero. The total optical conductivity in Matsubara space is 


(T"(iw„) = 2g ( 


2 ,C? + 1 


- du) 


(fip 




^ ,lV(p)(iV(p) + l) ■ * 


1 




dip) f 1 _ 

2p3 \2p + iUn 2p — iUr 


= (^- 


(C31) 


where ctq* is optical conductivity of a free scalar Held (Eq. A8 with m = 0). 
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